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Abstract 

Diffractive photoproduction of dijets was measured witli the ZEUS detector at 
the ep colhder HERA using an integrated luminosity of 77.2 pb~^. The measure- 
ments were made in the kinematic range 

Q2 < lGeV^ 0.20 < y < 0.85 and 
xip < 0.025, where is the photon virtuality, y is the inelasticity and xip is the 
fraction of the proton momentum taken by the diffractive exchange. The two jets 
with the highest transverse energy, -E;^*, were required to satisfy -E;^* > 7.5 and 
6.5 GeV, respectively, and to lie in the pseudorapidity range —1.5 < rf^^ < 1.5. 
Differential cross sections were compared to perturbative QCD calculations using 
available parameterisations of diffractive parton distributions of the proton. 
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1 Introduction 



In diffractive electron-proton scattering, the proton loses a small fraction of its energy 
and either emerges from the scattering intact, ep —>■ eXp, or dissociates into a low-mass 
state N, ep — > eXN. A large rapidity gap (LRG) separates the hadronic system X with 
invariant mass Mx and the final-state proton p. 

In the framework of Regge phenomenology [1], diffractive interactions are ascribed to the 
exchange of a trajectory with vacuum quantum numbers, the Pomeron trajectory. In 
quantum chromodynamics (QCD), the diffractive factorisation theorem [2,3] states that 
the diffractive cross section in deep inleastic scattering (DIS) can be expressed as the 
convolution of universal partonic cross sections and a specific type of parton distribution 
functions (PDF), the diffractive PDF (dPDF). Diffractive PDFs are interpreted as the 
number density of partons conditional on the observation of a diffracted proton in the final 
state. The dPDFs [4-7] have been determined from the HERA inclusive measurements of 
the diffractive structure function [4, 5] , defined in analogy with the proton structure 
function F2, and were used for input to calculations of hard diffractive processes at HERA, 
Tevatron and LHC [8-13]. 

Diffractive collisions producing a state X with a dijet system are a particularly interesting 
component of diffractive ep interactions. The transverse energies of the jets provide a hard 
scale, ensuring the applicability of perturbative QCD at the small photon virtualities con- 
sidered here. In photoproduction, at leading order (LO) of QCD, two types of processes 
contribute to dijet photoproduction, namely direct and resolved photon processes. In 
direct photon processes, the exchanged photon participates as a point-like particle, inter- 
acting with a gluon from the incoming proton (photon-gluon fusion. Fig. [1^). Thus, these 
processes are directly sensitive to the gluon content of the diffractive exchange. In resolved 
photon processes, the photon behaves as a source of partons, one of which interacts with 
a parton from the diffractive exchange (Fig. [lb). For resolved photon processes, which 
resemble hadron-hadron interactions, QCD factorisation is not expected to hold [3,14]. 
Further interactions between partons from the photon and the proton may fill the rapidity 
gap, leading to a breakdown of hard-scattering factorisation and causing a suppression of 
the diffractive photoproduction cross section. Such a mechanism was proposed to explain 
the suppression of the measured cross sections for hard diffractive hadron-hadron scatter- 
ing at the Tevatron with respect to expectations based on dPDFs obtained at HERA [15]. 
For the diffractive resolved photoproduction, an eikonal model [16] predicts a cross-section 
suppression by about a factor of three. In the framework of another model [17], assum- 
ing that diffractive collisions reflect the absorption of an incident particle- wave, it has 
been argued that the strong factorisation breaking observed in diffractive hadron-induced 
processes should not be seen in photon-induced processes. 
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This analysis presents measurements of the diffractive phot opro duct ion of dijets using the 
ZEUS detector at HERA. A 30-fold increase in luminosity was achieved compared to the 
previous ZEUS analysis [18]. This, in combination with the addition of a new forwarcEl 
detector, allows measurements to be made in a wider kinematic range. Differential cross 
sections based on these measurements are compared to next-to-leading-order (NLO) QCD 
predictions at the hadron level. The comparisons are also made separately for subsamples 
enriched with direct and resolved photopro duct ion. A similar study has been recently 
published by the HI Collaboration [9]. 



2 Experimental set-up 

This measurement is based on the data taken with the ZEUS detector at the ep collider 
HERA in 1999-2000 when electrons or positrons of 27.5 GeV were collided with protons 
of 920 GeV. The sample used for this study corresponds to an integrated luminosity 
£ = 77.2 pb~^ (12.1 pb~^ and 65.1 pb~^ for the e~p and e~^p samples, respectively)!^. A 
detailed description of the ZEUS detector can be found elsewhere [19]. A brief outline of 
the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [20], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle 
region 15° < 6 < 164°. The transverse-momentum resolution for full-length tracks is 
(^{Pt)/pt = 0.0058j9r © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAG). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test beam conditions, are a{E)/E = 0.18/v^ for electrons 
and a{E)/E = 0.35//E for hadrons, with E in GeV. 

In 1998, the forward plug calorimeter (FPC) [22] was installed in the 20 x 20 cm^ beam 
hole of the FCAL, with a small hole of radius 3.15 cm in the centre to accommodate the 
beam pipe. The FPC increased the forward calorimetric coverage by about one unit in 
pseudorapidity to r/ < 5. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 

^ From now on, the word "electron" will be used as a generic term for both electron and positron. 
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The backing calorimeter (BAG) consists of proportional tube chambers placed in the gap 
of the iron yoke. In the present analysis it was used in conjunction with the CTD and 
the GAL to identify cosmic muons that traversed the yoke. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'yp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [23], a 
lead-scintillator calorimeter placed in the HERA tunnel a.t Z = —107 m. 



3 Kinematics and reconstruction of variables 

Diffractive photoproduction in ep scattering (Fig. [1]), 

e{e)+p{p)^e{e')+X{X)+p{p'), 

is described in terms of the four-momenta of the incoming and scattered electrons, e 
and e', of the incoming and scattered protons, p and p', and of the hadronic system, X. 
The following kinematic variables are defined: the photon virtuality, = —q"^, where 
q = e — e', the squared photon-proton centre-of-mass energy, = {p + qY, and the 
fraction of the electron energy transferred to the proton in its rest frame (inelasticity), 

p ■ q W"^ 

y= 7^ — • 

p ■ e 2p ■ e 

The reaction can be considered to proceed through the interaction of the virtual photon 
with the diffractive exchange (Pomeron, IP). This process is described by the invariant 
mass, Mx, of the hadronic system X and the fraction of the proton momentum carried 
by the diffractive exchange 

{p -p')-q 

Xip = . 

p-q 

In the present data, the state X contains a dijet system as the result of a hard scattering 
process. The partons from the resolved photon and the diffractive exchange participating 
in the interaction have fractional momenta given by 

p ■ u 
p ■ q 

where u is the four-momentum of the parton in the resolved photon, and 

V ■ q 
[p -p')-q 
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where v is the four-momentum of the parton in the diffractive exchange. 

Energy flow objects (EFOs) were reconstructed from CAL clusters and CTD tracks and 
combine the CTD and CAL information to optimise the resolution of the reconstructed 
kinematic variables [24]. The EFOs were additionally corrected for energy loss due to 
inactive material in front of the CAL [25]. 

The mass Mx of the hadronic system X was reconstructed as 

Mx = VE. (E - pz)h ■ EhiE + Pz)h , 

where the index h runs over all EFOs. The quantities E and pz indicate the energy and 
the longitudinal momentum of the EFOs, respectively. 

The inelasticity, y, was estimated from the EFOs according to the Jacquet-Blondel method 
[26] as 

yM = EhiE-pz)j2E, , 

where Eg is the initial electron energy. For events with an electron candidate, the inelas- 
ticity was also determined from the scattered electron, yg. 

The longitudinal momentum fraction transferred from the proton to the diffractive ex- 
change, xip, was reconstructed as 

xiP = EhiE + Pz)k/2Ep , 

where Ep is the initial proton energy. 

The jets were reconstructed from the EFOs by using the algorithm [27] in the lon- 
gitudinally inclusive mode [28] in the laboratory frame. The variables i?^*^'^ and rf'^^^''^ 
characterise the two jets with highest transverse energy, Et, with E}^^^ > For the 

variables x^ and zjp, which are not measurable directly, the observable estimators x"^^ [29] 
and 0?p^ were reconstructed as 



obs ^ Z^jetl,2 ^ 

2yE, 

^obs _ Z^jetl,2 ^ 

where the sums run over the two highest E^ jets. 

In direct-photon processes, at LO in QCD, x^ is equal to one, whereas resolved-photon 
processes appear at < 1. A direct-enriched region was defined by x"^*^ > 0.75 and a 
resolved-enriched region by x°}''^ < 0.75. 
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4 Event selection 



A three-level trigger system was used to select events online [19,30]. Events with a large 
energy deposit in the calorimeter, neglecting the three inner rings of cells around the 
beampipe in the FCAL, were selected at the first-level trigger. Additional cuts were 
applied at the second-level trigger to reject beam-gas interactions and other non-ep back- 
ground events. At the third level, the measured transverse energy, excluding the first inner 
ring of the FCAL, was required to be greater than 11 GeV. Jets were not preselected at 
any trigger level. 

Well-reconstructed events were selected by applying the following quality cuts. The events 
were required to have at least three well-measured tracks of transverse momentum pt > 
0.2 GeV originating from the same vertex. The longitudinal position of the vertex, Z^txi 
had to be in the range —35 < Z^tx < 30 cm. 

Photoproduction events were selected as follows. Events with a scattered electron candi- 
date having an inelasticity of ye < 0.7 were assumed to be DIS events and removed. In 
addition, 0.20 < i/jb < 0.85 was required. The cut on and the upper cut on yjB reduced 
the remaining background from DIS events and also restricted the range of the virtuality 
of the exchanged photon to < 1 GeV^ with a median value of 10^'^ GeV^. The lower 
cut on ?/jB removed proton-beam gas events which deposit energy in the FCAL near the 
beam pipe. 

Events with at least two jets were selected by requiring a transverse jet energy above 
> 7.5 (6.5) GeV. Both jets were required to be in the pseudorapidity range —1.5 < 
i^jeti,2 ^ measured in the laboratory frame. 

Diffractive events were selected by requiring the presence of a LRG between the scattered 
proton and the rest of the hadronic final state. Since the proton was not measured, 
the requirement of a LRG was implemented by a cut on the total energy in the FPC, 
Eppc < 1.0 GeV, and by demanding r^^ax < 2.8. Here ?7max is defined as the pseudorapidity 
of the most forward EFO with an energy above 400 MeV in the CAL. This selection ensures 
at least a two-unit rapidity gap in the hadronic system, suppressing background from non- 
diffractive and proton-dissociative processes. In addition, a cut xjp < 0.025 was applied 
to enhance the Pomeron-exchange contribution [31]. 

Finally, cosmic-ray events originating from muons that traversed the detector near the 
interaction point were removed. Details can be found elsewhere [32,33]. A total of 6990 
events remained after all selection cuts. 
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5 Monte Carlo simulations 



Monte Carlo (MC) simulations were used to determine acceptances and resolution effects 
at detector level and to extract the hadronisation corrections for the NLO predictions, 
i.e. ratios of event yields at hadron level to those at parton level. 

The MC generator Rapgap [34] was used to simulate dijet processes in diffractive pho- 
toproduction at the Born level. Electroweak radiative effects were simulated by using 
Rapgap in conjunction with Heracles 4.6 [35]. 

The electron-proton interactions at small were modelled with both direct and resolved 
photon processes (Fig. [1]). Events were generated assuming that diffractive processes 
proceed via the emission of a particle-like Pomeron from the proton followed by the 
interaction of the virtual photon with the Pomeron. Although this factorised approach 
has no justification in QCD, it gives a fair description of the data. The diffractive PDFs, 
as determined by the HI Collaboration (HI LO fit 2) [36] for the Pomeron contribution, 
were used. For resolved photon processes, the photon PDFs GRV-G-HO [37] were chosen. 

In the simulation chain, the process of QCD radiation is followed by hadronisation. This 
was simulated by interfacing Rapgap to a parton-shower model as implemented in Meps 
[38] and to a hadronisation model based on string fragmentation [39,40] as implemented 
in Jetset [41]. 

The generated MC events were passed through the standard simulation of the ZEUS 
detector, based on the Geant program [42], and a trigger-simulation package [30]. The 
simulated events were reconstructed and selected in the same way as the data. 

Since the MC events generated with Rapgap did not adequately describe the dis- 
tribution of the data, they were reweighted to the measured distribution separately for 
^obs ^ g ^obs ^ 0.75. The relative fractions of direct photon and resolved photon 

processes were determined from a fit to the data. Resolved processes account for about 
one third of the total event sample. 

Event distributions are compared with the reweighted Rapgap MC distributions for 
the kinematic variables Mx, E}^^'^ and r/j''*^'^ in Fig. [2j The MC distributions were 
normalised to the data yielding a reasonable overall description of the data. 

The hadronisation corrections were calculated with the Rapgap MC sample after reweight- 
ing its parton level z'^^ distribution to each of the NLO predictions described in Sec- 
tion 18.21 In addition, hadronisation corrections were also calculated with a MC sample 
generated with POMWIG [43], a modification of the Herwig MC program [44] based on 
a cluster fragmentation model [45,46]. Since only direct photon interactions can be simu- 
lated with PoMWiG, the comparison to Rapgap was restricted to the range x°}'^ > 0.75. 
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The bin-by-bin differences between the corrections obtained with the two programs give 
an indication of the systematic uncertainties due to the hadronisation corrections [33]. 

The MC generator Pythia [47] was used to model the non-diffractive photoproduction 
of two jets. Events were generated using the CTEQ5L [48] (GRV-G-HO) parametrisation 
of the proton (photon) PDFs and processed through the same simulation and selection 
chain as the data. 

6 Background 

Background from proton-dissociative events, with a low-mass proton-dissociative system 
escaping down the beam pipe, was estimated to be (16 ± 4)% [32] by fitting to the FPC 
energy distribution, without the Eppc cut, a mixture of Rapgap and Epsoft MC [49]. 
This value was also obtained from hard diffractive production of open charm [50]. It was 
assumed that this estimate is independent of all kinematic variables studied here. The 
measured cross sections were scaled down accordingly. 

Background from non-diffractive dijet photoproduction, as estimated with the Pythia 
MC, was found to be less than 5% throughout the whole kinematic range, and was ne- 
glected. 

7 Systematic uncertainties 

Systematic uncertainties on the measured cross sections were estimated as described be- 
low: 

• the trigger efficiency was estimated for both data and Monte Carlo events using an 
independent trigger branch. The efficiency was above 98% for the entire kinematic 
range. The Monte Carlo simulation agrees with the data within ±1% [33] and the 
uncertainty was neglected; 

• the transverse jet-energy scale was varied by ±3%, the typical uncertainty in this 
range [18]. This variation resulted in an uncertainty of less than ±5%; 

• the FPC energy cut was varied by ±0.5 GeV, resulting in an uncertainty less than 
±1% in most bins and not more than ±2%; 

• changing the energy threshold of the EFOs, which is used to calculate r^max, by 
±100 MeV led to an uncertainty typically less than ±1% and not more than ±2% 
in any bin; 
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• the ?7max values of data and Monte Carlo events were shifted relative to each other by 
±0.1, the typical t] resolution. This led to the largest observed uncertainties which 
were typically below ±6.5% and up to ±14% for low z'^p^ and large xjp and Mx', 

• the lower yjQ cut was varied within its resolution (0.04); the resulting uncertainties 
were typically less than ±1% and not more than ±3%. When the higher yjB cut was 
lowered, the measured cross sections changed typically by < 1% and not more than 
±4%; 

• varying the 77-'''* cuts within its resolution (0.1) gave an uncertainty which is mostly 
below ±1% and not more than ±4%; 

• the xip cut was varied within its resolution (0.0025); the resulting uncertainties were 
typically less than ±2%, increasing to ±6% in the highest Mx and rf'^^^ bins. 

The systematic uncertainties not associated with the jet-energy scale were added in 
quadrature to the statistical uncertainty and are shown as error bars of the measured 
cross sections in Figs [3] to O The uncertainty due to the energy scale is shown separately 
as a shaded band in each of the figures. Overall normalisation uncertainties of ±2.2% from 
the luminosity determination and of ±4% from subtraction of the dissociative background 
were not included. 



8 Results 

8.1 Cross sections 

Single-differential cross sections were measured in the kinematic region < 1 GeV^, 
0.20 < y < 0.85, xjp < 0.025, eJ.^"^^^ > 7.5 (6.5) GeV and -1.5 < 7]^"^^'^ < 1.5, and 
were determined function of y, Mx, xip, zp, E^^'\ r/j^" and xf\ The estimated 
contribution of proton-dissociative background of 16% was subtracted in all bins. 

The cross sections are shown in Figs [3] and H] and listed in Tables [IH3l The cross section 
dependence on x'^'^, shown in Fig. HI indicates that direct-enriched (x"''** > 0.75) processes 
dominate diffractive dijet photoproduction in the kinematic range of this measurement. 

Single-differential cross sections were also determined separately for direct-photon en- 
riched and resolved-photon enriched processes. They are shown in Figs [5] and [6l respec- 
tively, and listed in Tables |4] and [5l The two sets of distributions differ in shape. Typically, 
resolved events are characterised by larger diffractive masses Mx] this in turn reflects into 
the observed xjp behaviour. Slight differences are observed in the Zjp'^ distributions with 
the most prominent feature being the rise of the direct-enriched component when z'^}? 
approaches one. 
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8.2 Comparison to the NLO QCD calculations 

NLO predictions for diffractive photoproduction of dijets were calculated at parton level 
with a program by M. Klasen and G. Kramer [51]. The calculations were performed with 
a fixed-flavour number of Nf = 4 and A4 = 330 MeV, chosen to match the value of the 
running as in the region of four active flavours. Three sets of dPDFs were used: the 
ZEUS LPS fit, determined from an NLO QCD fit to inclusive diffraction and diffractive 
charm-production data [4], and the HI 2006 fits A and B, obtained from fits to inclusive 
diffraction data [5] . The Regge-inspired parameters set for the NLO calculations were the 
same as used to obtain the dPDFs. The t-slope used in the Pomeron flux was 5 GeV^. 
For comparison with data, the NLO calculations obtained with the HI dPDFs were scaled 
down by a factojf] of 0.87 [5]. The contribution of subleading Regge trajectories as imple- 
mented in the HI fits was included. For the resolved photon, the 7-PDF parametrisations 
GRV [37] and AFG04 [52] were used. 

The NLO QCD predictions were obtained setting the renormalisation and factorisation 
scales to = jj^F = jJ' = -E^j^*^. The theoretical uncertainties were estimated by varying 
the scales simultaneously between (0.5 ■ E}^^^) and (2 ■ E^^^^) [51]. Changing the number 
of active fiavours to Nf = 5 in the NLO calculations leads to an increase of the expected 
cross section for x"^^ > 0.75 by less than 10%, and to a negligible effect elsewhere. The 
uncertainties of the dPDFs and the Pomeron flux, constraining directly the normalisation, 
were not included. The predicted cross sections were transformed to the hadron level using 
the hadronisation corrections calculated with Rapgap as described in Section [51 The 
uncertainties of the hadronisation corrections are not included in the error calculations 
for the cross sections. 

The data are compared with NLO QCD predictions at hadron level for the full x?^^ range 
in Figs [3] and m The hadronisation corrections applied to the NLO predictions at parton 
level are shown in the lower part of each plot and the values are given in Tables [TH3l The 
asymmetric theoretical uncertainties, estimated as described, were determined for the 
ZEUS LPS flt; those for the other NLO predictions are similar. The data are reasonably 
well described in shape. However, they lie systematically below all the predictions. Most 
of the suppression originates from the lower Ei^^^ region. 

Figure Hl^b) shows the ratio of the data and the NLO predictions using the ZEUS LPS 
flt. The ratio is below one, consistent with a suppression factor of about 0.7 independent 
of x^. Also shown is the ratio expected if the calculated resolved-photon cross section is 
suppressed by a factor of 0.34 [16]. No additional suppression factor for resolved-enriched 

^ The HI measurements used to derive the HI dPDFs include low-mass proton-dissociative processes 
which leads to an overestimate of the photon-diffractive cross section by a factor of (l.tsig gg) as 
obtained from MC simulations [5]. 
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data is observed. The suppression factor depends on the dPDFs and ranges between 
about 0.6 (HI 2006 fit A) and about 0.9 (HI 2006 fit B). Within the large uncertainties 
of NLO calculations, the data are compatible with no suppression, as expected in [17]. 

Differential cross sections for the direct-enriched and resolved-enriched samples are com- 
pared with NLO predictions at hadron level in Figs \5\ and El Again the hadronisation 
corrections are shown in the lower part of each plot and the values are given in Tables H] 
andO For direct-enriched data, the hadronisation corrections are shown for both Rapgap 
and POMWIG. The differences are taken as an estimate of the uncertainties as described 
in Section [51 The data lie systematically below the NLO calculations. Also, contrary to 
NLO expectations, the cross section as a function of zjp^ for the direct-enriched sample 
rises steadily with increasing Zjp^. 

Compared to NLO calculations obtained with the program of Frixione and Ridolfi [53], the 
HI Collaboration observed a suppression factor of about 0.5 in both resolved-enriched and 
direct-enriched cross sections of diffractive dijet photoproduction [9] . The measurements 
of ZEUS and HI cover different kinematic regions in Et and x/pjj In particular, the HI 
measurements extend to lower values than in the ZEUS analysis. In ZEUS, the largest 
discrepancy between the measured and predicted values of the cross section is observed 
at the lowest Et values suggesting that the conclusion on factorisation breaking depends 
on the probed scale. 



9 Conclusions 

Cross sections for diffractive photoproduction of dijets were measured with the ZEUS 
detector at HERA using an integrated luminosity of 77.2 pb~^. The measurements were 
performed in the kinematic region Q2 < 1 GeV^ 0.20 <y < 0.85 and xip < 0.025. The 
two jets with highest transverse energy were required to have Ei^^^^'^^ > 7.5 (6.5) GeV and 
-1.5 < 7/j^*i'2 < 1.5. 

The measured differential cross sections are compared to NLO QCD predictions based on 
available parameterisations of diffractive PDFs. The comparisons were made for the full 
data sample as well as for the subsamples enriched with resolved photon {x"^^ < 0.75) 
and direct photon (x"*^^ > 0.75) processes. The NLO calculations tend to overestimate 
the measured cross sections of both the resolved-enriched and the direct-enriched data 
sample. However, within the large uncertainties of the NLO calculations the data are 
compatible with QCD factorisation. 



^ It was checked that both programs for calculating NLO predictions give consistent results. 
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y bin 


da/dy 


f^stat 


^syst 


^ES 


Chad 




(pb) 


(pb) 


(pb) 


(pb) 




0.20, 0.33 


137.9 


5.3 


+ 15.2 
-7.0 


+4.0 
-2.9 


1.16 


0.33, 0.46 


198.4 


6.3 


+18.9 
-5.0 


+12.1 
-4.1 


1.09 


0.46, 0.59 


218.3 


6.7 


+16.5 
-10.4 


+19.4 
-12.8 


1.05 


0.59, 0.72 


196.5 


6.2 


+19.1 
-8.2 


+13.0 
-5.7 


1.05 


0.72, 0.85 


203.6 


6.4 


+17.1 
-6.6 


+16.6 
-10.4 


1.04 


Mx bin 


da/dMx 


^stat 




5es 


Chad 


(GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




15.0, 20.0 


2.11 


0.11 


+0.11 
>-0.01 


+0.13 
>-0.01 


1.37 


20.0, 25.0 


6.35 


0.17 


+0.24 
-0.13 


+0.26 
-0.23 


1.16 


25.0, 30.0 


6.39 


0.17 


+0.68 
-0.15 


+0.51 
-0.26 


1.04 


30.0, 35.0 


5.41 


0.17 


+0.57 
-0.32 


+0.36 
-0.21 


1.01 


35.0, 40.0 


3.14 


0.14 


+0.62 
-0.19 


+0.32 
-0.13 


1.01 


40.0, 45.0 


1.21 


0.09 


+0.23 
-0.08 


+0.21 
-0.07 


0.97 


Xip bin 


da/dxjp 


^stat 


"^syst 


<^ES 


Chad 




(nb) 


(nb) 


(nb) 


(nb) 




0.000, 0.005 


0.40 


0.06 


+0.01 
-0.04 


+0.04 
>-0.01 


1.46 


0.005, 0.010 


3.94 


0.14 


+0.18 
-0.10 


+0.14 
-0.21 


1.21 


0.010, 0.015 


6.28 


0.17 


+0.36 
-0.10 


+0.53 
-0.20 


1.10 


0.015, 0.020 


7.00 


0.19 


+0.84 
-0.31 


+0.47 
-0.27 


1.02 


0.020, 0.025 


7.21 


0.21 


+1.13 
-0.50 


+0.62 
-0.18 


1.03 


zi^^ bin 


da/dzp 


f^stat 


^syst 


^ES 


Chad 




(pb) 


(pb) 


(pb) 


(pb) 




0.2, 0.4 


86.4 


5.0 


+12.4 
-9.0 


+9.9 
-3.6 


0.88 


0.4, 0.6 


145.7 


4.9 


+15.7 
-8.9 


+9.3 
-8.2 


0.92 


0.6, 0.8 


192.9 


4.9 


+9.6 
-8.0 


+11.3 
-10.7 


1.11 


0.8, 1.0 


190.2 


4.2 


+11.2 
-0.2 


+10.8 
-5.2 


1.49 



Table 1: Differential cross sections for the diffractive photoproduction of dijets as a 
function of y, Mx, xjp and zp listed with statistical (Sstat) md systematic fSsystj un- 
certainties and energy scale (6es) uncertainties; the last column shows the hadronisation 
corrections (C^ad) applied to the NLO QCD predictions. 
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E}^^^ bin 


da / dE}^^^ 


^stat 


'^syst 


5es 


Chad 


(GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




7.5, 9.5 


44.0 


0.8 


+3.9 
-1.6 


+0.4 
-1.0 


1.02 


9.5, 11.5 


13.7 


0.4 


+1.3 
-0.5 


+3.1 
-0.7 


1.22 


11.5, 13.5 


3.5 


0.2 


+0.3 
-0.2 


+0.4 
-0.3 


1.22 


13.5, 15.5 


0.8 


0.1 


+0.1 
-0.0 


+0.1 
-0.1 


1.35 




da/dr]^^^^ 


^stat 


•^syst 


<^ES 


Chad 




(pb) 


(pb) 


(pb) 


(pb) 




-1.5, -1.0 


14.2 


0.9 


+1.5 
-0.6 


+1.0 
-1.0 


0.70 


-1.0, -0.5 


55.6 


1.8 


+3.1 
-2.4 


+3.2 
-1.9 


0.94 


-0.5, 0.0 


72.9 


2.0 


+5.0 
-3.6 


+4.7 
-3.7 


1.09 


0.0, 0.5 


63.6 


1.8 


+6.3 
-2.2 


+4.5 
-2.6 


1.16 


0.5, 1.0 


34.9 


1.3 


+4.9 
-0.8 


+2.3 
-0.2 


1.28 


1.0, 1.5 


8.0 


0.6 


+1.1 
-0.3 


+0.9 
-0.3 


1.38 



Table 2: Differential cross sections for the diffractive photoproduction of dijets as a 
function of E^l^^^ and rf^^^ listed with statistical fJstaty' CLiT'd systematic fSsysty' uncertainties 
and energy scale (6-es) uncertainties ; the last column shows the hadronisation corrections 
(Ctiad) applied to the NLO QCD predictions. 



a;°^^ bin 


da/dx^l"' 
(pb) 


^stat 
(pb) 


^syst 
(pb) 


(pb) 


Chad 


0.250, 0.375 
0.375, 0.500 


28.5 
52.7 


2.4 
3.2 


+1.4 
-3.0 
+2.7 
-4.1 


+0.5 
-1.1 
+4.2 
-1.9 


1.23 
1.04 


0.500, 0.625 


78.1 


3.7 


+3.2 
-4.2 


+1.9 
-4.5 


1.01 


0.625, 0.750 


114.3 


4.5 


+6.8 
-7.4 


+5.9 
-7.3 


1.18 


0.750, 1.000 


356.5 


6.2 


+29.1 
-14.1 


+23.3 
-14.4 


1.07 



Table 3: Differential cross sections for the diffractive photoproduction of dijets 
as a function of x'^^ listed with statistical f^statj one? systematic (5syst) uncertain- 
ties and energy scale (5-e,s) uncertainties; the last column shows the hadronisation 
corrections applied to the NLO QCD predictions. 
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y bin 


da/dy 


^stat 


'^syst 


5es 


Chad 




(pb) 


(pb) 


(pb) 


(pb) 




0.20, 0.33 


123.2 


5.1 


+12.0 
-6.7 


+2.3 
-3.8 


1.10 


0.33, 0.46 


152.6 


5.7 


+13.0 
-4.5 


+10.2 
-4.8 


1.07 


0.46, 0.59 


151.8 


5.7 


+11.0 
-7.3 


+13.6 
-8.5 


1.07 


0.59, 0.72 


125.2 


5.0 


+13.4 
-5.3 


+7.9 
-1.8 


1.05 


0.72, 0.85 


135.3 


5.4 


+11.3 
-4.9 


+10.6 
-9.8 


1.09 


Mx bin 


da/dMx 


^stat 






Chad 


(GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




15.0, 20.0 


2.08 


0.11 


+0.10 
-0.01 


+0.13 
>-0.01 


1.37 


20.0, 25.0 


5.27 


0.16 


+0.30 
-0.08 


+0.22 
-0.13 


1.11 


25.0, 30.0 


4.88 


0.16 


+0.38 
-0.17 


+0.35 
-0.35 


1.05 


30.0, 35.0 


3.26 


0.14 


+0.36 
-0.25 


+0.23 
-0.14 


1.03 


35.0, 40.0 


1.70 


0.11 


+0.36 
-0.16 


+0.19 
-0.06 


1.05 


40.0, 45.0 


0.58 


0.07 


+0.12 
-0.07 


+0.12 
-0.07 


1.05 


xjp bin 


da/dxjp 
(nb) 


f^stat 


'^syst 


^ES 


Chad 


0.000, 0.005 


0.40 


0.06 


+0.01 
-0.05 


+0.03 
>-0.01 


1.48 


0.005, 0.010 


3.33 


0.12 


+0.18 
—0.05 


+0.14 
—0.15 


1.22 


0.010, 0.015 


4.84 


0.15 


+0.27 
-0.09 


+0.34 
-0.16 


1.10 


0.015, 0.020 


4.86 


0.16 


+0.56 
-0.26 


+0.34 
-0.22 


1.06 


0.020, 0.025 


4.43 


0.18 


+0.68 
-0.40 


+0.34 
-0.21 


1.03 


zp bin 


da/dzp 


^stat 


^syst 


^ES 


Chad 




(pb) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




0.2, 0.4 


72.0 


4.7 


+9.1 
-8.1 


+7.4 
-3.3 


0.89 


0.4, 0.6 


105.3 


4.2 


+12.3 
-5.6 


+6.5 
-7.1 


1.02 


0.6, 0.8 


120.6 


3.8 


+9.5 
-3.1 


+8.0 
-4.8 


1.17 


0.8, 1.0 


144.0 


3.8 


+4.4 
-2.4 


+7.3 
-5.2 


1.59 



Table 4: Differential cross sections for the diffractive photoproduction of dijets for 
^obs > Q listed with statistical f^statj o,i^d systematic ((^systj uncertainties and en- 
ergy scale (6es) uncertainties; the last column shows the hadronisation corrections 
applied to the NLO QCD predictions. 
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1 • 

y bill 


da /ay 


r 


r 

'-'syst 


^ES 


'-'had 




(pb) 


(pb) 


(pb) 


(pb) 




0.20, 0.33 


14.7 


1.5 


+2.5 

-i!o 


+0.7 

-o!i 


1.93 


0.33, 0.46 


44.9 


2.7 


+5.6 
-2.7 


+0.1 
-1.7 


1.19 


0.46, 0.59 


66.5 


3.5 


+11.9 
-6.9 


+3.9 
-5^9 


1.01 


0.59, 0.72 


71.0 


3.6 


+22.4 
-4.4 


+3.6 
-5.0 


1.08 


0.72, 0.85 


69.5 


3.5 


+18.6 

-4.5 


+4.3 
-2.8 


1.09 


Mx bm 


da /dMx 






r 


(--had 


(GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




15.0, 20.0 


0.03 


0.01 


+0.01 
-0.01 


<+0.01 
>-0.01 


1.72 


20.0, 25.0 


1.09 


0.07 


<+0.01 
-0.15 


+0.04 

-o!ii 


1.54 


25.0, 30.0 


1.55 


0.07 


+0.27 
-0.07 


+0.06 
-0.04 


1.07 


30.0, 35.0 


2.14 


0.10 


+0.59 
-0.18 


+0.07 
-0.14 


1.01 


o r n /inn 
00. U, 4U.U 


i.4o 


u.uy 


+0.50 
-0.14 


+0.07 
-0.12 


111 
1.11 


40.0, 45.0 


0.64 


0.06 


+0.26 
-0.06 


+0.07 
-0.04 


1.08 


Xip bin 


da 1 dxjp 


Ostat 


"syst 


^ES 


'-'had 




/ 1 \ 

(nb) 


/ 1 \ 

(nb) 


/ 1 \ 

(nb) 


nb) 




0.005, 0.010 


0.60 


0.05 


<+0.01 
-0.07 


<+0.01 
-0.07 


1.16 


0.010, 0.015 


1.44 


n no 

0.08 


+0.20 
-0.09 


+0.13 
-0.08 


1.13 


0.015, 0.020 


2.19 


0.10 


+0.46 
-0.13 


+0.04 
-0.14 


1.08 


0.020, 0.025 


2.76 


0.12 


+0.99 

-0.28 


+0.15 
-0.11 


1.13 


z'Ip'' bin 


da/dzp 


^stat 


"^syst 


^ES 


Chad 




(pb) 


(pb) 


(pb) 


(pb) 




0.2, 0.4 


14.4 


1.8 


+8.5 
-1.5 


+1.8 
-0.9 


0.96 


0.4, 0.6 


40.4 


2.6 


+22.3 
-3.7 


+2.5 
-1.4 


0.95 


0.6, 0.8 


72.3 


3.1 


+20.7 
-5.6 


+3.8 

-5.5 


1.19 


0.8, 1.0 


46.2 


1.8 


+0.8 

-11.3 


+0.8 
-2.7 


1.26 



Table 5: Differential cross sections for the diffractive photoproduction of dijets for 
< 0.75 listed with statistical (6stat) o-nd systematic (S^yst) uncertainties and en- 
ergy scale f^Esj uncertainties; the last column shows the hadronisation corrections 
applied to the NLO QCD predictions. 
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(a) (b) 

Figure 1: Leading-order diagrams for (a) direct and (b) resolved processes in 
diffractive photoproduction of dijets at HERA. The variables shown in the plots are 
described in the text. 
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Figure 2: Comparison of the data (dots) with the Rapgap MC (solid line) 
normalised to the data as a function of (a) y, (b) Mx, (c) E^^'^ , (d) 7]^""'^ (e) E^'*^ 
and (f) rj^^^'^ after all cuts but the one on the plotted variable. 
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Figure 3: Single- differential cross sections (dots) as a function of (a) y, (b) 
Mx, (c) xip, (d) zjp^ , (e) E}^^^ and (f) rf'^^^ compared with NLO Q CD predictions, 
corrected for hadronisation, using the dPDFs from the ZEUS LPS fit (solid line), 
the HI 2006 fit A (dashed line) and the HI 2006 fit B (dotted line) and the GRV 
'j-PDFs. The inner error bars of the dots show the statistical uncertainty, the outer 
error bars show the statistical and systematic uncertainties (see Section^ added in 
quadrature. The dark shaded band indicates the jet energy scale uncertainty. The 
light shaded band shows the theoretical uncertainty due to the variation of the scale 
when using the ZEUS LPS fit. Underneath each plot the hadronisation corrections 
applied to the NLO prediction at parton level are shown. 
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Figure 4: (a) Single- differential cross section as a function of x"^*^ compared 
with NLO QCD predictions, corrected for hadronisation, using the dPDFs from the 
ZEUS LPS fit (solid line), the HI 2006 fit A (dashed line) and the HI 2006 fit B 
(dotted line) and the GRV 'y-dPDF. The prediction with HI 2006 fit A is also shown 
using the AFG parametrisation of the •y—PDFs (dashed-dotted line). Other details 
are the same as in the caption of Fig. O (b) Ratio of data and NLO predictions 
using the ZEUS LPS fit and GRV. The histogram indicates the expectation with the 
predicted resolved photon component scaled down by a factor of 0.34. The shaded 
and hatched bands show the theoretical uncertainty. Underneath the hadronisation 
corrections applied to the NLO prediction at parton level are shown. 
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Figure 5: Single- differential cross sections as a function of (a) y, (b) Mx, (c) xjp 
and (d) z°p^ for direct-photon- enriched dijet photoproduction (x"^^ > 0.75) com- 
pared with NLO QCD predictions, corrected for hadronisation, using the dPDFs 
from the ZEUS LPS fit (solid line), the HI 2006 fit A (dashed line) and the HI fit 
B (dotted line) and the GRV '-/-PDFs. Underneath each plot hadronisation correc- 
tions are shown which were obtained with Rapgap (upper histogram) and POMWIG 
(lower histogram) , respectively. The shaded bands indicate the differences. The cor- 
rections from Rapgap were applied to obtain the NLO predictions shown above. 
Further details are the same as in the caption of Fig. 
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Figure 6: Single- differential cross sections as a function of (a) y, (b) Mx, (c) xjp 
and (d) z'^p^ for resolved-photon- enriched dijet photoproduction (x"^'^ < 0.75) com- 
pared with NLO QCD predictions, corrected for hadronisation, using the dPDFs 
from the ZEUS LPS fit (solid line), the HI 2006 fit A (dashed line) and the HI 
2006 fit B (dotted line) and the GRV •y-PDFs. Underneath each plot the hadronisa- 
tion corrections applied to the NLO predictions at parton level are shown. Further 
details are the same as in the caption of Fig. 0. 
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